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Electromagnetic form factors of proton and neutron, obtained from a new fit of data, are presented. 
The proton form factors are obtained from a simultaneous fit to the ratio [i p Ge p /Gm v determined 
from polarization transfer measurements and to ep elastic cross section data. Phenomenological 
two-photon exchange corrections are taken into account. The present fit for proton was performed 
in the kinematical region Q 2 £ (0,6) GeV 2 . Both for protons and neutrons we use the latest 

, available data. For all form factors the uncertainties and correlations of form factor parameters are 

■ investigated with the \ 2 method. 
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1— 1 1. INTRODUCTION 

The nucleon electromagnetic form factors are fundamental quantities, of great theoretical and experimental impor- 
tance. The issue of their determination has been revisited in recent years, thanks to the results of several experiments 
at Bates, MAMI, JLab, which put under question previous analyses based on less precise data and urged the ne- 
cessity for a new parameterization and a new analysis of the form factors themselves (for a review, see, for example 
Refs. 

A precise knowledge of the electromagnetic form factors of the nucleon is important for the determination of the 
Q\ ' axial nucleon form factor in charged current (CC) quasielastic neutrino-nucleon scattering Q and strange form factors 
C*~) . of the nucleon in neutral current (NC) elastic neutrino-nucleon scattering. For example, NC vector form factors which 
I characterize elastic NC scattering are given by the following expressions Q : 

. G^ C,P '"' = ±- {Ge p — Geu} — 2 sin 2 6wGEp(n) ~ ~2*~* Es 

£,NC,p(n) _ _|__ {Q Mp _ Gmu} — 2 Sin 2 OwG^ip^n) — -Gms ■ 

• »"H , In the above the dominant terms are the electric (Ge) and magnetic (Gm) form factors of the nucleon. Their precise 
knowledge is essential in order to determine the small strange form factors of the nucleon Ge s , Gms- Hence, it is 
obvious that not only a good knowledge of the electromagnetic form factors is required, but also the present level of 
" ' their uncertainty. 

In this paper we performed new fits of the nucleon electromagnetic form factors. The proton ones are extracted 
from: i) elastic ep cross section data, ii) polarization data, providing the ii p Ge/Gm ratio (/i p being the magnetic 
moment of the proton). The neutron form factors are extracted from electron-nucleus (typically deuterium and 3 He) 
scattering processes. The latest experimental data are used. 

The proton form factors determined from the measurements of polarization transfer in elastic electron-proton scat- 
tering (first appearing between '99 and '02) were in a significant disagreement with respect to the ones obtained from 
elastic ep scattering data via the customary Rosenbluth separation. The main suggestion to solve this inconsistency 
was to account for two photon exchange (TPE) diagrams jy, 0, B @| , which should affect the cross section to a greater 
extent than the polarization data. 

This disagreement became even more evident after the new JLab data on ep scattering cross sections |10| appeared. 
Hence, as already pointed out by several authors (see e.g. [ill]), a reliable global fit must include the TPE correction; 
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we will explicitly show the effect of TPE on the goodness of the fit (GoF). For a recent review devoted to TPE 
correction see Ref. [12fl . 

There exist several parameterizations of the nucleon form factors which have been considered in the literature 
JT3JT3 JT1 EE EE E3, Hi, H HE [2l|, [22)]. Among these, the older ones have a purely empirical Q 2 dependence ES 

EJEIEJ 



G Ep (Q 2 ), G Mp (Q 2 )/v P jf , Q = VQ 2 - 

1 + c ^ 



i=l 



The specific form of the parameterization may depend on the Q 2 region. For instance, in Ref. [19( the low-Q 2 data 
were analyzed with form factors given by a continued fraction parameterization: 



G E P (Q 2 ) , G m p (Q 2 ) I ' i^p 



1 



hQ 2 



1 



1 + 



The newest empirical form factors are constrained to have a proper physical behavior at low-Q 2 as well as at high-Q 2 . 
One of the examples is the Kelly's parameterization EH, which will be employed in our analysis (see next Section). 
The form factors depend on powers of the invariant Q 2 , and for large Q 2 the form factors behave like l/Q 4 . In 
Ref. [20( Kelly's parameterization is additionally constrained to satisfy duality hypothesis and the low-Q 2 behavior is 
described as in Ref. EH- 

The electric neutron form factor, usually, is separately treated and described with a smaller number of parame- 
ters [i| (see also [3). 

It is also worth mentioning those parameterizations obtained on the basis of the vector meson dominance model. In 
particular the parameterization proposed by Lomon [22j seems especially suited to successfully describe the neutron 
form factor data. 

In the present paper we aim to provide reliable fits of both proton and neutron e.m. form factors, by employing a 
relatively small number of parameters; moreover one of the major merits of this work is the analysis of errors on the 
parameters of the fit, which allows one to estimate the present uncertainty on our knowledge of the electromagnetic 
form factors. 

The paper is organized as follows: in Section 2 we consider the proton form factors, by analyzing both the polar- 
ization data (Section 12. ip and the cross section data (Section 12. 2|) . Section 3 is devoted to the neutron form factors, 
separately considering the electric and the magnetic form factors. Finally Section 4 presents a discussion of our results 
in comparison with previous analyses and the conclusions. 

2. PROTON FORM FACTORS 

In Section 12.11 we consider the recent polarization transfer and asymmetry measurements data which give an 
information on the ratio of the electric and magnetic proton form factors. Then, in Section [2. 2[ we present the results 
of the combined fit of the polarization and cross section data. 

2.1. Fit of polarization data 

In this section we consider the direct determination of the ratio of the electric and magnetic proton form factors 

which has been obtained with the measurement of the polarization of the recoil proton and with asymmetry mea- 
surements. Here Q 2 = — q 2 , q being the four- momentum transfer. In the one-photon approximation q is the four- 
momentum of the virtual photon. 
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FIG. 1: Linear fit ((9} of the recoil polarization and asymmetry measurements of the ratio /-l p Gep/Gm p . The shadowed area 
denotes the 3<r C.L. region of the fit. 



The recoil polarization technique (see Ref. [3j) has been employed in several ep experiments for a direct measurement 
of the ratio 1Z(Q 2 ). In the laboratory frame it is given by: 



K(Q 2 ) 



P t E + E ' 
2M 



tan 



(2) 



where Pi and P t are the longitudinal and transverse components of the polarization of the recoil proton, M is the 
proton mass, E and E' are the initial and final electron energies, and 9 is the electron scattering angle. The latter is 
related to Q 2 according to: 

Q 2 = AEE'sm 2 (- 



(3) 

The ratio 1Z(Q 2 ) has been also determined from the measurement of the asymmetry in elastic ep scattering with 
both polarized beam and target: we have (see Ref. 



= -2/i p yr(i + r) tan - 



0\ ft sin 0* cos 0* +fi p ^T[l + (1 + t) tan 2 (§)] 



cos 6* 



H 2 + n pT /e 



(4) 



where a + and cr_ are the cross sections for positive and negative electron helicities, respectively, 9* and <f>* are the polar 
and azimuthal angles of the target polarization relative to the three-momentum transfer vector q and the scattering 
plane (in the laboratory frame) , 



and 



Q 2 

AM 2 ' 



l + 2(l + r)tan^ - 



(5) 



(6) 



is the virtual photon polarization. 

We consider the recoil polarization and asymmetry data published in Refs. m m m m m nn m m, m, 



35|, l3fj|. These data are plotted in Fig.[I] together with their error bars, which include the statistical and systematic 
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FIG. 2: Projection on the a — /3 parameter space of the contours delimiting the allowed regions in the 2-dimensional space of 
TPE correction parameters, with la (dashed lines) and 2a (solid lines) C.L. The contours are computed for fit I (left figure) 
and fit II right figure. The crosses indicate the projections of the best-fit point, Eq. (|28[1 for fit I and Eq. (|3f I) for fit II. 



uncertainties added in quadrature. The Q 2 range of the data goes from 0.15 to 5.6 GeV 2 . As one can see from Fig. [IJ 
all data are well described by a linear function in Q 2 : 

K(Q 2 )=c + Cl Q 2 . (7) 

We fitted the data points with this linear function, by minimizing the least-squares function 



N »* (TZ(Q 2 ) - Kf p ) 2 



where iV rat = 65 is the total number of recoil polarization and asymmetry data points and 7^ xp is the value of the 
ratio at the squared-momentum transfer Q 2 , with corresponding uncertainty A7?° xp . 
We found the following best fit values of the parameters: 

c = 1.022 ±0.005, ci = -0.130 ±0.005, (9) 

with la uncertainties computed from the covariance matrix 1 (they are given by the square-roots of the diagonal 
elements of the covariance matrix). 
The corresponding minimum \ 2 being: 

(X? at )min/NDF = 58.89/63, (10) 



where NDF = A ra t is the number of degrees of freedom; the goodness of the fit (see Ref. (38|) is 62%. 

The solid line in Fig. [1] corresponds to the best-fit values of the parameters in Eq. while the shadowed area 
denotes 3cr C.L. region of the fit. One can see that the linear fit has small uncertainties, especially for Q 2 < 3 GeV 2 , 
where there are many data points. 

The best-fit values of the parameters Co and c\ in Eq. (El ar e close to those obtained by Arrington in Ref. [la ]. 
c{j f = 1.0324 and c\ l = -0.135. Let us notice that in Refill it was assumed that for Q 2 < 0.24 GeV 2 the form 
factor ratio is equal to one. 



2.2. Fit of polarization and cross section data 

The values of the proton form factors have been extracted from the data of many elastic ep scattering experiments 
using the Rosenbluth method. In the one-photon approximation, the differential cross section in the laboratory frame 



1 Detailed numbers of the covariance matrices relative to these and to the following fit parameters can be found in Ref. |37l| . 
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FIG. 3: The fits of Gm p / h p Gd and Ge p /Gd form factors. The solid lines denote fit I. The shadowed areas represent the 2<r 
(bright) and 3u C.L. (dark) allowed regions. Kelly's fit [l8l ] is shown by the dotted line. 



for unpolarized ep elastic scattering reads (in the same notation used in the previous subsection): 

da 1 -v ,„„t,/„o t „o \ / 1 

e 

ctm being the Mott's differential cross section 



^^^-m(^Q 2 )( G | p + ^L p ) — . (ID 



^ Q ] = [d^e) M - 2£3 s in>/2) ■ (12) 
The Rosenbluth separation is then obtained by considering the reduced differential cross section 

a K {E, Q 2 ) = e (1 + r) ^fj,^ = r G 2 (Q 2 ) + e G|„(Q 2 ) . (13) 



A linear fit of the reduced differential cross section at fixed Q and different values of e gives the value of tG m p (Q ) 
from the intercept (e = 0) and the value of G 2 Ep {Q 2 ) from the slope. Notice, however, that the measurement of 
G 2 Ep {Q 2 ) with the Rosenbluth method has large uncertainties, because the contribution of G 2 Ep {Q 2 ) to the reduced 
differential cross section in Eq. (fT3")l is suppressed for large values of Q 2 (r > e) while for small values of Q 2 we have 
G e p — G m p /^p — g m p / 7 - 8 - 

In our analysis, in the first fit, later called fit I, we assume that Ge p is related to Gmp by the linear relation of 
Eq. 0, which is favored by the direct measurement of TZ(Q 2 ) in polarization experiments, as discussed in Section l!01 
For the proton magnetic form factors we adopt the parameterization proposed by Kelly [l8[ : 



Gm p (Q z 



1 _l \ " „M 
1 + 2^ a P,k T 
k=l 

n+2 

M 



(14) 



J2K 

k=l 



which guarantees the asymptotic behavior Gmp(Q 2 ) °c Q 4 at high Q 2 [3§|. We shall employ the parameterization 
of Eq. (fT4"| with 4 parameters (n = 1): 



G Mp (Q 2 ) _ l + « 



T 



II 1 _1_ h M n- _1_ h M J 2 _L ^ M ^3 

Mp 1 + O pl T + b p2 r +0 p3 T 



(15) 
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We will see that this choice turns out to be quite satisfactory for the description of the data. Moreover a relatively 
small number of parameters allows a better control of the errors. 

We have also performed a fit with both the magnetic and electric proton form factors parameterized by the expression 
(fT5)) . This fit will be called fit II in the following. In this case the electric form factor reads: 

G Ep (Q 2 ) = — Te— ~ h E P '\ , rR 3 ( 16 ) 

1 + b p,l T + b P,2 T + b p,3 T 



In our analysis we consider similar sets of cross section data as the ones employed by Arrington in Ref. 17j , namely 
the data from Ref. [4(|-[63]. Some of the data were taken from the JLab data base [63| ; we include also data from 
Ref. [64| • Additionally, we considered the latest data of JLab experiment [l(j in which the cross section was measured 
with the smallest errors, up to-date. 

We fitted the ep cross section data by minimizing the least-squares function 



where M cs = 24 or M cs = 28 2 are the numbers of data sets, Nf s is the number of points in the ith data set, rii and A 



a, 



are the corresponding overall normalization and uncertainty, a i • is the jth differential cross section point in the ith 
data set, with electron energyi^^ and four-momentum transfer Q 2 j, a~(Ei ; j, Q 2 •) is the corresponding differential cross 
section computed with Eq. pip . The uncertainty Ac^- of c°j P includes the statistical and uncorrelated systematic 
uncertainties added in quadrature. 

We perform a simultaneous fit of the polarization and cross section data by minimizing the sum of the least-square 
functions, Eqs. © and p7| : 

X 2 = X?at + Xcs • (18) 

For the fit I the range of Q 2 taken into account is 

0.1 GcV 2 < Q 2 < 6GeV 2 , (19) 

which corresponds to the interval of Q 2 values where polarization transfer data are available. 
For the fit II we extend the range of Q 2 down to Q 2 ~ 0: 

< Q 2 < 6 GeV 2 . (20) 

Notice that a precise knowledge of the form factors in the low Q 2 region is of special interest for neutrino-nucleon 
(-nucleus) scattering processes. In both cases the upper limit, Q 2 = 6 GeV 2 , is determined by the polarization data. 
We do not consider higher Q 2 points. 

As already stressed in the literature [llj, the inclusion of the most precise data of Ref. [l(| definitely indicates 
the need of corrections to formula (fTTj) . Beyond the classical radiative corrections [65], to get agreement with the 
polarization data, one needs to consider also the two photon exchange (TPE) corrections, which can be written as an 
additive term to the reduced cross section: 

o~r — > o R + Stpe- (21) 

The calculation of 5tpe is difficult and model dependent: however, one can use general properties to derive a 
phcnomenological expression of the TPE term. The scattering amplitude for electron-nuclcon interaction must satisfy 
general symmetry properties, such as crossing symmetry and C-invariance [66j. They can be used to constrain Stpe- 



Following Ref. [67| we adopt a TPE correction given by a function F(Q 2 , y) 

cn^<7 R + TF(Q 2 ,y) (22) 



2 For Q 2 G (0.1, 6) we have 24 independent data sets while for Q 2 S (0, 6) we have 28 independent data sets. 
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FIG. 4: Electric and magnetic proton form factors. The solid lines denote the best fit II. The shadowed areas represent the 
2a (bright) and 3<r C.L. (dark) allowed regions. Our previous fit I is shown with the dashed lines. 



where 



V = 



1 - e 
T+e 



(23) 



The function F(Q 2 ,y) must satisfy the relation F(Q 2 ,y) — —F(Q 2 , —y). The analytical properties of F(Q 2 ,y) allow 
one to express this function as a series of odd powers of y. Chen et al. [(57j truncated the expansion to the second 
term: 



F(Q 2 ,y) = aG 2 D (Q 2 )y + [3G 2 D {Q 2 )y\ 
a and (3 being fit parameter and Gd(Q 2 ) the usual dipole form factor: 

G D {Q 2 )=[l + ^j , with My = 0.71 GeV 2 . 



(24) 



(25) 



We consider both types of fit; for the fit I we obtained: 

xLn/NDF = 375.97/392, GoF = 71% , 
with the following values for the best fit parameters: 

af x = 1.53 ±0.01, bf } \ = 12.87 ±0.07, bf 2 = 29.16 ±0.25, = 41.40 ± 0.33 , c = 1.02 ±0.01, c x = -0.13 ±0.01. 

(27) 

The parameters of the TPE correction are: 

a = -0.39 ±0.09, (3 = -0.04 ± 0.09 . (28) 



(26) 



Notice that the values of c and C\ parameters are very similar to the ones in Eq. @, obtained by fitting the 
polarization transfer data alone. 

For the fit II the minimization procedure leads to: 



xLn/NDF = 450.95/468, GoF = 71% , 



(29) 
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FIG. 5: The ratio [i v Ge v IGm v obtained in the simultaneous fit to polarization measurements and cross section data, fit I is 
shown by the dashed line, the corresponding 3<r C.L. allowed region being denoted by dark shadowed area, fit II is shown by 
the solid line, the corresponding 3<r C.L. allowed region being denoted by bright shadowed area. 



with the following values for the fit parameters: 

bf %x = 12.31 ± 0.07 , b^ 2 = 25.57 ± 0.22 , b™ 3 = 30.61 ± 0.27, a£f x = 1.09 ± 0.01 , 

bp A = 11.12 ± 0.15 , = 15.16 ± 1.03 , = 21.25 ± 3.27 , a^ = -0.19 ± 0.06. (30) 

The parameters of the TPE correction are: 

a = -0.36 ±0.09, (3 = -0.08 ± 0.09. (31) 

We remark that from both fits we obtained comparable values of the TPE parameters (see Fig. [2j which illustrates 
the allowed regions in the (a,/?) parameter space with a given confidence level (C.L.)). In both cases the TPE 
correction turns out to be negative. Let us mention that the way we introduce the TPE corrections in our analysis 
also motivates the choice for the upper Q 2 limit: indeed, following the approach of Ref. [67|], the magnitude of TPE is 
fitted to the data and, in the elastic cross section, it can be comparable to the magnitude of Ge p - Hence the inclusion 
of the polarization data (which are less affected by TPE correction) allows a more precise determination of the TPE 
fit parameters, but restricts the Q 2 range to the one of the available polarization data. 

It is worth mentioning that by excluding the TPE correction (hence using for the cross section formula (fTTj) ) both 
fits worsen, particularly in the goodness of fit. For the fit I we obtain x^ in /NDF = 467.07/394 with GoF = 0.6%; 
similarly for the fit II we get x^m/NDF = 544.31/470 with GoF = 1%. We noticed that this result stems from the 
presence, in the analysis, of the very accurate JLab data [10|, without which GoF would increase to 45% and 47%, 
respectively. 

In addition to the above discussed form factors, we also checked a different parameterization, based on a two-poles 
formula for both the electric and magnetic proton form factors [j|: 

GeviQ 2 ) = i , a L +7^ ( 32 ) 
Gm p (Q 2 ) 





af 


1 


-af 


1 - 


f afQ 2 


1 + 


aiQ 2 




af 


1 


-af 


1 - 


f afQ 2 


f l4 


-*fQ 2 



(33) 



With respect to Kelly's parameterization this would offer the advantage of having a smaller number of parameters, 
in addition to the ones of the TPE correction. A new, global fit, can be obtained with Xmin/NDF = 1.10 but 
GoF = 0.06%, thus indicating that the former parameterization is preferable. 
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FIG. 6: Comparison of fits obtained with and without TPE correction (solid and dashed lines correspondingly). In the first 
and second panels Gm p /(^ p Gd) and Ge p /Gd are plotted, respectively. The shadowed areas denote the 3<r. allowed regions. 
In the bottom panel the ratio (|37[) is shown. Results are obtained for fit II. 

As it has been already mentioned in the introduction, one of our main tasks is to compute the form factor uncer- 
tainties as they can be extracted from the fit. This goal can be achieved by performing an accurate error analysis on 
the various fit parameters. 

We calculated the correlated uncertainties of the fit parameters and the related uncertainties of the form factors with 
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the standard least-squares method, which is appropriate and widely used for non-linear models 3 (see Ref. 38l.l68l.l6 



the allowed region in the space of N parameters with A confidence level (C.L.) is delimited by the contour defined by 

X 2 =xLn + A X 2 (N,X), (34) 

where Ax 2 (N, A) is the value for which a \ 2 variable with N degrees of freedom has a cumulative probability A. We 
consider 2a (95.45% C.L.) and 3cr (99.73% C.L.) uncertainties. 

Since we have 8 parameters in fit 1 and 10 parameters in fit II, the exploration of the parameter space in order to 
find the contours defined by Eq. (f34|) cannot be done with the simplest grid method. Therefore, we used a Monte Carlo 
Markov Chain generator of random points, which allows to find the allowed parameter regions with good accuracy in 
a few hours of CPU time of a normal PC. 

It is interesting to notice, that the estimated values of the magnetic form factor parameters, a^, b p ! A , b^f 2 and b^f 3 
are strongly correlated. In particular, the estimates of a\ and 63 are almost linearly dependent. These parameters 
determine the asymptotic behavior of Gmp(Q 2 ), which turns out to be: 



M 

lim Q i G Mp {Q 2 ) = (4Af 2 ) 2 ~ (0.681°;°, GeV 2 ) 2 , (35) 



in fair agreement with the one given by the usual dipole form factor (|25|) (the above uncertainties are at 3a). 

In Fig. [2] (left panel) the error contours for a and (3 parameters are shown. Let us notice that solutions with (3 
positive but very small are possible, but in this case a should be negative and large in magnitude. Therefore the TPE 
correction are always negative. 

In Fig. [3] we show our best fits (fit I) for the magnetic and electric proton form factors with their uncertainties. We 
compare with Kelly's fit [lij ]. One can see that Kelly's fits of the magnetic proton form factor lies within our 3a C.L. 
region in almost the whole Q 2 range under consideration. For the electric form factor the fits differ by more than 3a 
in a relatively wide range of Q 2 . 

A similar error analysis is performed for fit II. Here, the number of form factor parameters is larger (4 parameters 
for each form factor). Similarly as above we show contour plot for the TPE correction parameters (Fig. [2 right panel). 

Analogously to the case of fit I the estimates of the parameters and b^ 3 are linearly dependent and their ratio 
is similar: 



lim Q i G Mp {Q 2 ) = (4M 2 ) 2 * (0.66±°;°2 GeV 2 ) 2 . (36) 

Q^oo b p3 

One could also derive from the form factor parameters the charge and magnetic root-mean-square radii for the 
proton, as given by the slope of the electric and magnetic form factors at Q 2 = 0. They turn out to be: J < r 2 Ep > = 



0.87 ± 0.01 fm, yj < r 2 Mp > = 0.86 ± 0.01, fm. These results are comparable with previous analysis in the literature, 
but slightly lower than the most recent and advanced estimates of Ref. 7(], 71 1 ; indeed the latter take into account 



Coulomb distortion, which is relevant at low Q 2 . For example Ref. [7l| provides yj < r 2 Ep > — 0.895 ± 0.018 fm. 

The present fit of the form factors is carried out in plane-wave approximation and low-Q 2 properties like charge 
and magnetic radii are not properly reproduced without accounting for radiative corrections to the Rosenbluth cross 
sections. 

In Fig. [4] we present Gmp and Ge p obtained in fit II, with the 2a and 3a C.L. error bands (represented by 
shadowed areas). Here the results from fit I are also plotted. Both fits lead to very similar magnetic form factors. On 
the contrary, there is a visible difference between the corresponding electric form factors: the Ge p obtained in fit I 
decreases faster then the one obtained in fit II. 

Fig. O shows the ratio i^ p Ge p /Gm p obtained with the fits I and II. The linear ratio fitted only to the polarization 
data is no longer shown, since it is very similar to the one obtained with fit I. The ratio uncertainties are larger for the 
fit II than for fit I due to the fact that the parameterization in fit II contains a larger number of degrees of freedom. 



In this case, the uncertainties of the parameters and their correlations estimated from the covariance matrix are quite approximate. This 
method was widely used in the past, when computer power was insufficient to perform more accurate evaluations, as the one presented 
in this paper. 
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FIG. 7: Fit of the electric neutron form form factor (solid line). The Galster-like parameterization (|38[) is considered, 
shadowed areas denote the 2a (bright) and 3a C.L. (dark) allowed regions. Kelly's fit pl| is shown with the dotted line. 



The 



Finally, given for granted that they are necessary, it is interesting to understand which is the quantitative impact 
of the TPE correction: they are expected to be relevant especially for the electric form factor. For this purpose we 
compare in Fig. [S] the proton form factors obtained with and without TPE correction - only fit II is considered. The 
magnetic proton form factor obtained without TPE is systematically shifted down by about 1.5%, and above Q 2 ~ 1 
GeV 2 it lies outside the 3cr C.L. region of the form factor obtained by including the TPE correction. The analogous 
effect on the electric form factor is shown in the middle panel of the same figure and appears to be less uniform than 
for the magnetic form factor: this can be better appreciated from the bottom panel of Fig. [6l where the ratio 



G Ep (without TPE)/G_e p (with TPE) 



(37) 



is plotted. One can see that the TPE correction substantially alters the Q 2 dependence of the electric form factor, in 
particular, for Q 2 > 2 GeV 2 , with an effect which grows up to the order of 10%. In any case the impact of the TPE 
correction turns out to be non-negligible for both form factors. 

3. NEUTRON FORM FACTORS 

The measurement of the neutron form factors is much more difficult than that of the proton form factors, since 
a target of free neutrons does not exist. The neutron form factors are extracted from measurements of electron- 
nucleus scattering, usually electron-deuteron or electron-helium scattering. Therefore, the data analysis is affected by 
uncertainties stemming from the nuclear theoretical model assumed to describe the target nucleus. Since these models 
have consistently improved with time, in our analysis we consider only relatively recent data. At variance with the 
proton case, we take from the literature directly the published values of neutron form factors "data" and apply our 
fitting procedure to them. 



3.1. Electric neutron form factor 



For the electric neutron form factor we adopt the Galster-like parameterization 



Geu(Q 2 ) 



At 
1 + Bt 



G D {Q 2 ). 



(38) 
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FIG. 8: Fit of the electric neutron form form factor (solid line). The BLAST-like parameterization (|42|l is shown. The shadowed 
areas represent the 2<r (bright) and 3a C.L. (dark) allowed regions. The dashed line denotes the Galster-like parameterization. 
The Lomon result (GKex02S) [22j is denoted by the dash-dotted line. 




CT [GeVl 



FIG. 9: Fit of magnetic neutron form factor (solid line), normalized to the dipole form factor. The shadowed areas denote 2<r 
(bright) and 3<r C.L. (dark) allowed regions. Kelly's fit [18| is shown with dotted line. 



with the dipole form factor of Eq. (f2"5|) . 

We consider the electric neutron form factors "data" which have been published in several papers. Some of the 
data have been obtained in asymmetry and recoil polarization measurements 72 , 

M H MM BH IzI M MM, 

83, 84, 85|. We consider also the reanalyzed electron-deuteron data and the newest BLAST measurements [87]]. 
Additionally, in order to have a proper slope of the electric form factor in the limit Q 2 — > 0, we impose to our fit the 
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additional constrain [88|: 

(r| n ) = -0.1148 ± 0.0035 fm 2 . (39) 

We considered a least-squares function similar to the one in Eq. l[5]). with the experimental statistical and system- 
atical uncertainties added in quadrature. With the values: 

A = 1.68 ± 0.05 , 5 = 3.63 ±0.39, (40) 

we obtained 

xLn/NDF- 25.82/37, (41) 

and the goodness of the fit turned out to be excellent: 91%. 

In Fig. we plot the best-fit value of Geu as a function of Q 2 together with the 2a and 3tr C.L. allowed regions. 
We plot also Kelly's fit 0. 

As an alternative to the most commonly used Galster-like parameterization, we considered a neutron electric form 
factor given by the sum of two dipole form factors: 

G En {Q 2 ) = j x. (42) 

(I + &1Q 2 ) 2 (1 + 6 2 Q 2 ) 2 

This parameterization is similar to the one considered in the latest BLAST data analysis [87] and for this reason we 
will call it BLAST-like parameterization. The fitting procedure for the above parameterization leads to: 

X 2 /NDF = 17.95/36, GoF = 99% (43) 

with the parameters: 

a = -0.10 ±0.02, bi = 2.83 ±0.37, b 2 = 0.43 ± 0.11. (44) 

In Fig. [5] the BLAST-like parameterization is compared to the data: the Galster-like parameterization and one of the 
recent Lomon parameterization [2l|, [22| (GKex02S) are also shown. Notice that the parameterization (|44|) raises faster 
with Q 2 than the Galster-like and the Lomon one, but the latter remain both within the BLAST parameterization 
uncertainties. 



3.2. Magnetic neutron form factor 

For the neutron magnetic form factor we adopted again the simplest form of Kelly's parameterizations, with n = 1: 

G M n(Q 2 ) _ 1 + <1T 



Mn l + ftElT+fi&^+ft&T*' (45) 

We considered 11 data sets, obtain ed fro m asymmet ry m easurements [§!|-[92j], [sH and cross section measurements in 
electron-deuterium scattering 94|- |l03j . where Ref. [l03| contains the latest JLab measurements. The fit to all these 
data sets leads, how ever, to a minimum % 2 /NDF = 2.05, not quite satisfactory. According to a remark of Kelly [70l |. 



the data from [10l| and [97[ were extracted using the same associated-particle technique for the neutron efficiency, 
a technique which appears to be in contradiction with the method used in other experiments. Therefore we omitted 
these two data sets. From our final analysis. After excluding the two above mentioned data sets, we obtained a fit 
over N = 56 points with 

XLJNDF = 52.79/52 = 1.01, (46) 
and GoF = 44%. The corresponding values for the parameters in the neutron magnetic form factor are: 

6f 1 = 21.30 ± 4.56, 6f 2 = 77 ± 31 b^ 3 = 238 ± 105, = 8.28 ± 3.89. (47) 
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FIG. 10: Magnetic and electric proton form factors (normalized to the dipole form factor). The comparison between fit I (solid 
lines) and Arrington el al. [ll[ (dashed lines) is shown. The shadowed areas denote the 3<r C.L. allowed region. 



We performed the error analysis over the four parameters above. Even if the fit of Gmu was done on a slightly 
different basis than the one of the proton, yet we observe strong correlations between estimated values of the param- 
eters, in full analogy with our findings for the proton form factor. In particular, the estimated values of parameters 
(ai and 63) which determine the asymptotic behavior of the form factor at large Q 2 are almost linearly dependent: 

M 

lim Q 4 G M „(Q 2 ) - (4M 2 ) 2 ~ (0.66 ±°"°} GeV 2 ) 2 (48) 

This value is very similar to the one obtained for the proton. Notice, however, that without the newest JLab data, 
instead of the value gSJ) we would get (0.58 ±° ° 4 5 GeV 2 ) 2 . In Fig. [Hour final fit of G Mn is different from Kelly's 
result 1 181 since it contains the newest JLab measurements. 



4. DISCUSSION AND CONCLUSIONS 



In this section we start by presenting further comparisons of the form factors resulting from our fits with the ones 
of previous data analyses. 



In Ref. [llj the first systematical global analyses of the cross section and polarization transfer data on the proton 



with the inclusion of the TPE correction was performed. That fit is valid up to Q 2 = 30 GeV 2 for the magnetic 
form factor and up to 6 GeV 2 for the electric form factor. In Fig. [TO] we display together our global fits and those of 
Ref. [lH: it clearly appears that, even thought different approaches for the TPE correction were employed, the global 
fits are very similar. 

We also compare our fits with the recent one of Bodck et al. 20] (BBBA07). This global fit is tailored to accurately 
describe the form factors at low Q 2 as well as in the intermediate region of Q 2 . These authors used the Kelly's 
parameterization with four parameters but each form factor was multiplied by some Legendre polynomial, which 
depends on several additional parameters, constrained to reproduce the low Q 2 behavior obtained in Ref. [l9(. 

The authors of Ref. [2(| plotted the form factors against the so-called Nachtman variable, which for the elastic 
scattering is defined as £ = 2/(1 + yl + 1/r). Therefore in order to make the comparison with their results we 
express our form factors in terms of the £ variable. Plots are shown in Fig. 1111 In the region of £ corresponding to 
the range of validity of our fits the predictions of the two parameterizations are very similar, however, our magnetic 
proton form factor is systematically higher (by several percent) then the one given by the BBBA07 parameterization. 
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FIG. 11: Proton and neutron electric and magnetic form factors (normalized to the dipole form factor). The fits of electric and 
magnetic form factors (fit II) as well as the fit of neutron magnetic form factor are denoted by solid lines. BBBA07 20] global 
fits are denoted by dashed lines. In the case of the electric neutron form factor, both parameterization, Galster and BLAST 
-like, are shown by the solid and dash-dotted lines, respectively. The form factors are plotted against the Nachtman variable. 
The dashed vertical bar represents the upper bound of the Q 2 range were our fits were performed. 



The difference is given by TPE correction which we considered in our fitting procedure (see Fig. [5]), while authors 
of Ref. did not discuss this effect. For higher £ values one can notice sizeable deviations for our form factors. 
However, we notice, that the £ variable compresses in a very short range the large Q 2 region. 

Finally we compare our estimates of uncertainties with those obtained by Arrington and Sick 19]: these authors did 
a serious attempt to compute the uncertainties of the nucleon form factors, which is a crucial information in the study 
of parity violating ep scattering. Their method to compute errors is explained in Ref. [l9| and differs from our, in 
particular in the treatment of the systematic uncertainties. In particular we usually obtain asymmetric uncertainties 
around the best fit value: hence in Fig. [12] we compare the errors of Arrington and Sick on the electric and magnetic 
proton form factors with our lower and upper bounds for the 3cr confidence level errors. For the magnetic proton form 
factor, within the 3er C.L. our results are consistent with the ones of Ref. [n|. For the electric proton form factor we 
notice some deviations between our results and the ones obtained by Arrington and Sick. 

In conclusion we have presented two fits of the proton and neutron electromagnetic form factors, using the best 
available data. The ep elastic cross sections were reproduced by including a simple but realistic parameterization of 
the two photon exchange correction. Alternative parameterizations with fewer parameters than the one employed 
here do not allow to obtain equally good fits. We show that the impact of the TPE correction on the magnetic and 
electric proton form factors is larger than the 3<r uncertainty of the fits (in a wide range of Q 2 ). In fit I we constrained 
the electric proton form factor by the ratio \i v G e p / 'G Mp extracted from recoil polarization and asymmetry data, fit II 
employs Kelly's parameterization with four parameters both for the electric and magnetic proton form factors. This 
fit is obtained with two additional parameters with respect to fit I, however we believe that it is more reliable than 
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Q 2 [GeV 2 ] Q 2 [GeV 2 ] 

FIG. 12: The uncertainties of electric and magnetic proton form factors. Our computations (fit II) are denoted by dotted 
(upper uncertainties) and dashed lines (lower uncertainties) for the 3<r C.L. errors. The solid lines correspond to the results 
obtained in Ref. [jg| |. 



the former, particularly in the low Q 2 region. We also performed a careful analysis of the uncertainties resulting on 
the parameters of the fit and hence on the form factors. It is worth stressing that only a few papers, among the many 
devoted to the nucleon electromagnetic form factors, include the TPE correction in the analysis. As a final remark 
we remind the reader that even small uncertainties in the magnetic form factors of the proton and neutron turn out 
to be important for a correct analysis of the neutrino-nucleon cross sections. 



The numerical results of our fits are available in the web site [37[ . 
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